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Abstract 

Changes in ocean angular momentum about the polar axis (M) are related to fluctuations in 
zonal currents (relative component M r ) and latitudinal shifts in mass (planetary component 
Mu)- Output from a 1° ocean model is used to calculate global M r , Mq, and M time series at 
5-day intervals for the period January 1985-April 1996. The annual cycle in M T , Mq, and M is 
larger than the semiannual cycle, and Mu amplitudes are nearly twice those of M r . Year-to-year 
modulation of the seasonal cycle is present, but interannual variability is weak. The spectrum of 
M is red (background slope between w _1 and cj~ 2 ) at subseasonal periods, implying a white or 
blue spectrum for the external torque on the ocean. Comparisons with previous studies indicate 
the importance of direct atmospheric forcing in inducing subseasonal M signals, relative to 
instabilities and other internal sources of rapid oceanic signals. Regional angular momentum 
estimates show that seasonal variability tends to be larger at low latitudes but there are many 
local maxima due to the spatial structure of zonal current and mass variability. At seasonal 
timescales, latitudes ~20°S-10°N contribute substantial variability to Mq, while signals in 
M r can be traced to Antarctic Circumpolar Current transports and associated circulation. 
Variability in M is found to be small when compared with similar time series for the atmosphere 
and the solid Earth, but ocean signals are significantly coherent with atmosphere-solid Earth 
residuals, implying a measurable oceanic impact on length-of-day variations. 
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1. Introduction 

As a rotating geophysical fluid, the oceans possess 
angular momentum M about the polar axis consist- 
ing of a planetary component Mo, due to the solid 
body rotation, and a relative component M r , due 
to their zonal circulation. Study of the axial ocean 
angular momentum (OAM) is important to the un- 
derstanding of the variable Earth’s rotation. In the 
absence of external torques, the planet conserves its 
angular momentum. Thus changes in OAM or in 
atmospheric angular momentum (AAM) should be 
mirrored in records of the length-of-day (LOD). The 
close link between AAM and LOD at seasonal and 
shorter timescales has long been established, but loss 
of coherence appears at periods of about a month and 
amplitudes drop below significant levels at subweekly 
periods (see reviews by Hyde and Dickey [1991] and 
Rosen [1993] and papers therein). The ocean is thus 
expected to contribute at the residual level at seasonal 
timescales and become more important at higher fre- 
quencies. 

Knowledge of the oceanic mass and zonal veloc- 
ity fields and their variability is required to assess 
OAM and its role in closing the planet’s angular mo- 
mentum budget. Partly because of the lack of data, 
the first OAM studies focused only on relative com- 
ponent associated with individual currents such as 
the Antarctic Circumpolar Current (ACC) [Munk and 
MacDonald , 1960; also Ponte and Rosen , 1994, here- 
after PR94, and references therein]. Modeling devel- 
opments over the last decade, however, have produced 
a number of simulations of the global ocean circula- 
tion and thus provided a new tool to estimate ax- 
ial OAM [e.g., Brosche et al, 1990; PR94; Bryan , 
1997, hereafter B97; Ponte , 1997; Segschneider and 
Sundermann, 1997; Johnson et al, 1999]. Marcus et 
al [1998] have demonstrated the measurable impact 
of OAM signals on LOD fluctuations at timescales 
from 1 year to 10 days. Other studies [Ponte, 1990; 
Holloway and Rhines , 1991; Straub , 1993; Cummins , 
1995] have also used the angular momentum frame- 
work to examine aspects of the ocean circulation re- 
lated to the role of stress and pressure torques at the 
atmosphere and solid Earth interfaces. 

Despite the recent revived interest in the subject, 
knowledge of OAM continues to be sketchy and pri- 
marily focused on the seasonal variability of M r . Most 
of what is known comes from simple (forward) model 
simulations withoqt any data constraints and is inher- 
ently uncertain, although the comparisons with AAM 


and LOD are encouraging [PR94; Marcus et al , 1998]. 
Calculations of M r from simulated zonal currents de- 
pend on the model friction parameterizations, bottom 
topography representation, and surface forcing, as has 
been noted by B97 (e.g., available estimates of annual 
cycle in M r vary by a factor of two). Calculations of 
Mo, based on values of bottom pressure diagnosed 
from model density and sea level fields, are consid- 
erably more uncertain because of unrealistic density 
trends present in most models and because of the im- 
pact of the Boussinesq approximation (and in some 
models the rigid lid assumption) on the sea level fields 
[PR94; B97; Ponte , 1999]. New modeling studies are 
thus expected to provide further insight on M r and 
Mq signals and their uncertainties. 

Here we use output from a run of the newly de- 
veloped model of Marshall et al. [1997] to extend 
previous studies of axial OAM. Ponte et al [1998] 
and Ponte and Stammer [1999] have used the same 
output to demonstrate the role of the equatorial com- 
ponents of OAM in the excitation of polar motion, at 
the same time attesting for the quality of the model 
fields. After describing model and data sets in sec- 
tion 2, global OAM quantities are treated in section 
3, together with comparisons with previous results 
from other models that allow further insight into the 
uncertainties of model-based OAM estimates. Sec- 
tion 4 discusses regional OAM quantities and their 
relation to global values, in particular those related 
to Mq whose regional characteristics are for the first 
time revealed here. The impact of OAM signals on 
LOD, discussed first by Marcus et al [1998] using 
a different model, is revisited in section 5 and for a 
much longer period. A summary and some general 
remarks conclude the paper. 

2. OAM and Ancillary Datasets 

Our definitions and calculation of OAM quantities 
follow closely PR94 and Ponte et al [1998]. A fluid 
particle of density p and unit volume, rotating about 
the Earth’s polar axis, has axial angular momentum 
m given by 

m ” mo -f m r = pQr 2 cos 2 </> + pr cos <f> u (1) 

where r and fl are the mean Earth radius and rota- 
tion rate, respectively, <j> denotes latitude, u is zonal 
velocity relative to solid body rotation, and mo and 
m r are the planetary and relative components of the 
particle’s angular momentum. Integrating (1) over 
the ocean volume and using the hydrostatic relation 
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p z = — gp, with p being pressure and g being the ac- 
celeration of gravity, yields for the total planetary and 
relative OAM components 

Mfi = J j Pb COS 3 <f> d(j> d\ = ^ xf (2) 

M r — — J J j u cos 2 (f> dcj) dX dp = Cfixaf (3) 

Here A is longitude, p& is the pressure at the ocean 
bottom, and we have set the surface pressure to 
zero. The variables xf and Xsf are, respectively, the 
pressure and velocity excitation functions derived by 
Barnes et al. [1983] and commonly used in comput- 
ing the effects of the geophysical fluids on LOD. xf 
and X3" are basically equivalent to Mq and M r but 
nondimensionalized by CSl, where C = 7.04 x 10 37 
kg m 2 is the polar moment of inertia of the mantle. 
The extra factor of 0.7 in xf accounts for the effects 
of solid Earth deformation under loading. Thus, only 
70% of the Mq signals are effectively involved in LOD 
excitation. We will use x values when analyzing OAM 
in the context of the planet’s momentum budget and 
M r , Mq values otherwise, but conversion between the 
two quantities is straightforward using (2) and (3). 

Model output used to calculate Mq and M r is the 
same as in Ponte et al [1998] and Ponte and Stammer 
[1999]. For these studies, the model of Marshall et al 
[1997] was run on a near-global domain (80°S-80°N) 
with l°xl° horizontal grid spacing and 20 vertical lev- 
els. Step-like topography was used. The model was 
started from the final state of a 3-year spin-up run and 
forced with twice-daily surface wind stresses and daily 
surface heat and freshwater fluxes from the National 
Centers for Environmental Prediction, in addition to 
surface restoring to climatological temperature and 
salinity fields. More details on the model output can 
be found in the papers cited above and also Stammer 
et al [1997]. Integrals in (2) and (3) are carried over 
the model ocean domain, with p& values estimated as 
in Ponte [1999]. The spatially-uniform sea level cor- 
rection of Greatbatch [1994] was applied to account 
for global ocean volume changes due to unmodeled 
steric effects. With this correction, the total oceanic 
mass in the model is constant in time. 

To study the role of OAM signals on the planet’s 
angular momentum budget, values of AAM or \ A 
based on the National Centers for Environmental Pre- 
diction/National Center for Atmospheric Research re- 
analysis, are obtained from the Sub-Bureau for Atmo- 
sphere of the International Earth Rotation Service. 


Estimates of both \ V and X P are used. Contribu- 
tions from winds up to 10 hPa levels are included in 
X v ] values of x P are based on the inverted barome- 
ter assumption [Salstein et al , 1993], consistent with 
the timescales analyzed here [Ponte, 1997]. Geodetic 
excitation functions \ G are computed from SPACE97 
LOD time series [e.g., Gross , 1996]. 

All angular momentum and x series are calculated 
for more than 11 years (January 1985- April 1996) 
and are provided as 5-day averages. Besides globally- 
integrated OAM time series, we also calculate regional 
integrals over boxes of 1° x 1° (the horizontal resolu- 
tion of the model), which will be denoted as Mq, M r . 
The sum of gridded values over all boxes yields the re- 
spective globally-integrated values. However, local in- 
tegrals are based only on monthly-averaged quantities 
and thus capture mostly the seasonal variability. The 
gridded Mq, M r values are used to assess the impor- 
tance of different basins and regions in contributing 
to the seasonal cycle in OAM. No torque or momen- 
tum flux terms are available from this run, however, 
and no attempt is made to examine the OAM budget, 
either regionally or globally. 

3. Global OAM Variability 

Time series of Mq, M r , and their sum M in 
Figure 1 display variability at seasonal and shorter 
timescales, with weak interannual signals and no ap- 
parent secular trends. The annual cycle is clear in all 
time series, but with visible year-to-year amplitude 
modulations (e.g. M strong in 1989, weak in 1994). 
Signals in Mq tend to be larger than those in M r and 
noticeably so at the annual period. The importance 
of Mq is consistent with Ekman and geostrophic dy- 
namics that imply the ocean mainly adjusts to season- 
ally varying stress and pressure torques by changing 
its planetary angular momentum [PR94]. The rela- 
tive importance of Mq is in stark contrast with what 
happens in the atmosphere, where M r dominates sea- 
sonal variability [e.g., Salstein et al, 1993] and nonlin- 
ear eddy processes are important [Peixoto and Oort, 
1992]. The M r and Mq terms in Figure 1 are posi- 
tively correlated, with a coefficient of 0.69. Because 
of the positive correlation, variability in M is larger 
than that in either Mq or M r , with peak-to-peak am- 
plitudes ~10 25 kg m 2 s -1 or ~10% of the range of vari- 
ability in AAM [e.g., Rosen, 1993]. The correlation 
between M r and Mq is expected for geostrophic dy- 
namics, in which zonal currents are related to merid- 
ional pressure gradients and thus to latitudinal distri- 
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bution of mass. 

Spectra of Mq, M r , and M in Figure 2 show that 
power in Mq is larger than that in M r over all fre- 
quencies (but less so as frequency increases), con- 
tributing most of the power in M. Nevertheless, all 
spectra exhibit similar shapes and a background slope 
between u~ l and u )~ 2 . Such frequency dependence 
for spectrum of M implies a white or blue spectrum 
for the net external torque acting on the ocean, which 
is a sum of the surface zonal wind stress torque and 
east-west pressure and stress torques on the bottom 
boundary [Ponte, 1990; B97]. As the spectrum of the 
surface stress torque on the oceans is predominantly 
red (not shown), ocean dynamics play a fundamental 
role in determining the character of OAM variability 
resulting from atmospheric forcing. 

For a qualitative measure of the low frequency vari- 
ability, the periodograms of Mq , M r , and M are also 
shown in Figure 2. There is a tendency for whiteness 
at interannual timescales with a peak at the annual 
period. The overall importance of Mq in explain- 
ing changes in Af at all frequencies is clear. Despite 
the short records and the year-to-year modulation 
noted in Figure 1, we perform a harmonic analysis 
on 11 complete years of data (1985-95) to character- 
ize annual and semiannual cycles. Respective ampli- 
tudes and phases are given in Table 1. Amplitudes of 
Mq are nearly double those of M r , for both annual 
and semiannual periods, but nearly in phase (partic- 
ularly for the former). Annual cycle is almost twice 
as large as the semiannual cycle. Relative phasing 
of the annual and semiannual cycles leads to a clear 
maximum in OAM values in the late boreal summer 
(Figure 1). 

The present estimate of the annual cycle in M r has 
an amplitude of about half that of PR94, with a max- 
imum occurring approximately 1 month later, but is 
much more similar, both in amplitude and phase, to 
those reported by B97 for most of his experiments, ex- 
cept the one with smooth topography, and also with 
those in Brosche et al [1990]. For the semiannual 
cycle, amplitude spread is much wider, with values in 
Table 1 being about 2 to 3 times larger than those 
of B97 and PR94, and 7 times larger than those in 
Brosche et al. [1990]. With the exception of PR94, 
semiannual maximum occurs always in late Febru- 
ary to late March (and 6 months later). Overall, 
M r phase estimates seem to be more stably deter- 
mined than amplitudes. 

Plausible reasons for the range in seasonal M r val- 
ues include differences in forcing fields (e.g., PR94 


and B97 use climatological monthly mean winds), 
eddy parameterizations and resolution (e.g., PR94 es- 
timates are based on 0.5° x 0.5° grid, with better re- 
solved eddies), and bottom topography. The degree of 
smoothness of the topography field was found by B97 
to have the largest influence on the annual cycle of 
M r , compared to changes in wind fields and viscosity 
parameters. Marcus et al. [1998] also speculate on 
the role of topography representation in improving 
OAM estimates. We note that our model has compa- 
rable vertical resolution (20 levels) to that of PR94, 
but uses less smooth topography. The comparison 
with B97 values is thus consistent with the idea that 
rougher, more realistic topography may yield smaller 
annual amplitudes in M r . 

A comparison of present and previous Mq (and 
thus M) estimates is not as straightforward, because 
of the different assumptions made in dealing with the 
seasonal steric sea level component. We restrict dis- 
cussion to values of B97, who uses the same method 
to correct for unmodeled steric effects. Compared to 
B97 values computed with ECMWF wind climatology 
for 1986-92 (see his Table 4), amplitudes in Table 1 
are similar but somewhat larger for both the annual 
and semiannual cycles, and annual maximum occurs 
1 month later and semiannual maxima occur a couple 
of weeks earlier. Note that the different upper surface 
formulations in the two models (rigid lid in B97, free 
surface in our case) appears to have a relatively small 
impact on the calculation of Mq. 

Besides differences at the seasonal time scale, the 
subseasonal OAM variability in PR94 is much weaker 
than that shown in Figure 1. Time series of PR94 are 
based on an eddy-resolving model forced by monthly 
mean winds. The differences with the non-eddy re- 
solving, synoptically-forced estimates in Figure 1 thus 
suggest that intrinsic oceanic instabilities and related 
eddies are not as important as atmospherically-driven 
variability in producing rapid OAM signals, and that 
knowledge of wind stresses is essential in studies of 
such OAM signals ( Ponte [1997] makes the case for 
the importance of surface pressure at subweekly peri- 
ods). This finding contrasts with the apparent dom- 
inant contribution of intrinsic eddy processes to sur- 
face variability in many oceanic regions [e.g., Stam- 
mer, 1998]. Globally-integrated OAM of course fa- 
vors barotropic signals coherent on the largest scales, 
and these may be essentially atmospherically-driven. 
But eddies generated by oceanic instabilities may also 
carry little net angular momentum signals (either im- 
plying w T eak interaction with or nearly compensat- 
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ing torques on the bottom boundary). More definite 
answers to these issues must await calculations with 
fully-resolved eddies (eddy field in the model used by 
PR94 is only marginally resolved and its strength un- 
derestimated). 

4. Regional OAM Seasonal Variability 

An analysis of regional OAM variability can shed 
light on the importance of given circulation features 
to global OAM. The importance of u and pt, sig- 
nals involved in M r and A/q, respectively, increases 
with proximity to the equator where the moment arm 
a cos (/> is largest, but variability in u and pb also differs 
substantially from region to region (see Figure 8b in 
Ponte and Stammer [1999]) due to differences in ocean 
dynamics, forcing, bathymetry, etc. Our discussion of 
regional OAM signals at seasonal timescale comple- 
ments analysis of M r signals by PR94 and provides a 
first description of Mq signals, setting the stage also 
for future studies of regional OAM budgets, when 
local torques and flux terms will be available. We 
consider local integrals over 10° x 10° boxes to focus 
on the large scale. 

Figure 3a shows the standard deviation of M r , and 
the fractional covariance, defined as 

((Mr ~ (M r ))(M r - (Mr))) 

{{Mr - {Mr))*) 

where angular brackets represent time averaging, is 
shown in Figure 3b to highlight relation between lo- 
cal and global quantities. By its definition, fractional 
covariances over all grid boxes sum up to unity. Vari- 
ability in M r is strongest in the tropics, particularly 
in the Pacific and Indian oceans, where zonal circula- 
tion is known to have a strong seasonal modulation, 
and consistent with results in PR94 and B97. Much of 
the tropical variability, and more generally the basin 
variability, is, however, expected to recirculate at dif- 
ferent latitudes, as oceanic flows are to zeroth order 
horizontally nondivergent. The recirculation can be 
inferred from the alternating bands of negative and 
positive fractional covariances, strikingly aligned in 
the zonal direction especially in the Pacific, Indian, 
and Southern oceans. 

Major contributions to seasonal variability in M r are 
likely to come from the Southern Ocean where un- 
blocked zonal flow through Drake Passage is possible. 
The positive covariances at ~60°-20°S latitudes add 
up to nearly one and indicate the importance of in- 
cluding flows up to subtropical southern latitudes to 


capture full seasonal variability in M r (similar to find- 
ing by PR94). Knowing the ACC transport T through 
the Drake Passage, and assuming perfectly nondiv er- 
gent flow, one can estimate M r as 2n r 2 p co s 2 (f> T 
[Munk and MacDonald, 1960; PR94], where cos 2 <j> 
is based on some representative latitude. Figure 4 
shows such an estimate in comparison to the glob- 
ally integrated M r . The similarity between the two 
curves is clear at the seasonal timescale. (At shorter 
time scales, the assumption of horizontally nondiver- 
gent flows becomes less valid, but analysis not shown 
reveals significant coherence between Drake Passage 
transport and M r at all frequencies. The correlation 
coefficient between the two series is 0.66.) 

Variability in Mq, examined in Figure 5a, is weaker 
and not as spatially variant as that in M r . As before, 
the Pacific and Indian Oceans show larger signals than 
the Atlantic. Despite a tendency for weaker variabil- 
ity at high latitudes, there are local maxima in the 
Southern Ocean (west of South America and south- 
west of Australia) and in the western North Pacific. 
These regions were first noted for their large subsea- 
sonalp& fluctuations [Fukumori et al , 1998], but show 
similar enhancement at seasonal periods [Ponte, 1999; 
Ponte and Stammer , 1999]. Spatial variability of 
the wind stress curl or topographic features leading 
to a much reduced effective (3 may yield a locally en- 
hanced response [Ponte, 1999]. Variability in coastal, 
shallow seas (e.g., Sea of Japan, South China Sea, 
Arafura Sea between Australia and New Guinea) can 
be surprisingly large. These constricted regions are 
not optimally resolved in the model, but results seem 
to agree with the output of other models run at much 
higher resolution (cf. Plate 1 in Ponte [1999]). Large 
Pb (and thus Mq) signals in shallow areas can be ex- 
pected given that wind stress forcing terms scale in- 
versely with ocean depth. The large variability in 
equatorial regions may also be related to the strong 
monsoonal forcing, as well as the different dynamics 
of the equatorial wave guide. 

The relation between Mq and Mq is quite differ- 
ent than for M r . With few exceptions, mainly at high 
latitudes, positive covariances dominate in Figure 5b, 
indicating adding contributions to Mq over most re- 
gions. Covariances are, however, small (maximum be- 
low 0.02) in comparison with the case of M r in Figure 
3b. The contributions from the Atlantic are weak- 
est compared to the other oceans. Tropical latitudes 
(~20°S-10°N) contribute most strongly to seasonal 
Mq signals. The characteristics of Mq reflect the 
weak but very large scale nature of seasonal anomalies 
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in pb [Ponte, 1999]. 

With M r signals being stronger than those in Mq , 
the pattern of variability in M (Figure 6a) follows 
in general that of M r , with maxima in tropical In- 
dian and Pacific Oceans. These maxima coincide with 
maximum positive covariances in Figure 6b. Compar- 
ing Figures 3b, 5b, and 6b, M has a weaker banded 
structure than seen in Figure 3b because of the posi- 
tive contributions over most regions from Mq. Thus 
the structure of added covariances in longitude in 
Figure 6b resembles that in Figure 3b, but values 
are generally positive and most important in two 
bands at ~20°-40°S and ~10°-20°N. In the South- 
ern Ocean, Mq variability compensates somewhat for 
that in M r and decreases the importance of the latter 
at ~40°-50°S. 

5. OAM, AAM and LOD 

The impact of OAM signals in Figure 1 on the 
planet’s axial angular momentum budget is assessed 
by comparing them to observed AAM and LOD vari- 
ability. The motivation for this analysis is to check 
the degree to which our OAM estimates represent 
measurable signals in LOD. Attempting a “best” 
angular momentum budget would involve thorough 
analysis of all the time series involved, including AAM 
and LOD, and is beyond our intended scope. 

Figure 7 compares the oceanic atmospheric 
X3 , and geodetic x$ time series, and spectra and co- 
herence results are displayed in Figure 8. (Note that, 
in calculating x° > effects of Mq are reduced by 30%; 
see discussion in section 2.) As noted in the Intro- 
duction, time series of \z and X3 are closely linked 
and almost indistinguishable to the eye. Variability 
in residuals X 3 -X 3 ( or Xz~ A for short) is at around 
10% of that in either series, and most conspicuous at 
subseasonal periods. Spectra of X 3 and X3 have very 
similar shapes and power levels at ail periods resolved, 
consistent with a dominant relation between AAM 
and LOD. Power levels in x^ are, however, generally 
higher, indicating missing LOD excitation. (Largest 
discrepancies occur near the fortnightly band, where 
X3 7 spectrum shows some enhancement not matched 
by X3 spectrum; incomplete removal of solid Earth 
tide in the Xz series may be responsible for this fea- 
ture.) 

Figures 7 and 8 clearly show the weak variability in 
Xz compared to xf or X3 4 ■ Variance in Xz is only 6% 
of that in X3 , and also small (31%) compared to that 
in the residual x$ ~ A ' Spectrum of Xz shows weaker 


power than that in Xz ~ A spectrum at all frequency 
bands. Variability in OAM is, however, significantly 
coherent with the residual series Xz~ A most pe- 
riods (Figure 8), demonstrating the presence of cor- 
related signals in both series. Furthermore, adding 
OAM to AAM signals increases the coherence with 
LOD in general and most clearly at periods of 13-18 
days, providing also for slight improvements in the 
phase estimates (i.e., values closer to zero). 

Our results agree qualitatively with the findings of 
Marcus et al [1998] and confirm the observable role of 
OAM signals in the planet’s axial angular momentum 
budget. However, the amount of atmosphere-solid 
Earth residual variance explained by our estimated 
OAM signals is small (only 10%) compared to their 
best results (42% in the case of the Miami Isopycnal 
Coordinate Ocean Model run). Besides the different 
ocean models and periods analyzed, there are also two 
important differences in the way AAM series are com- 
puted. First, Marcus et al AAM values include wind 
contributions from the uppermost levels of the atmo- 
sphere (10-0-3 hPa), which seem to reduce substan- 
tially the seasonal x^ A residuals and may effectively 
improve agreement with OAM signals at seasonal pe- 
riod. Second, Marcus et al do not consider Mq con- 
tributions to AAM. We find that by similarly neglect- 
ing atmospheric Mq signals, Xz~ A variance explained 
by OAM signals rises to ~20%, However, in this case, 
variance in Xz ~ A also goes U P by more than 30% and 
high frequency coherence between Xz and Xz is S en_ 
erally lower than in Figure 8. (In fact, residual rms 
magnitude in Figure 7 is smaller than that reported 
by Marcus et al (52.7 ps vs. 60.5 ps using LOD 
units), and coherence amplitudes in Figure 8 are no- 
ticeably higher compared to their Figure 3b. Thus, 
it seems that including Mq signals in AAM provides 
better agreement with LOD.) For these reasons, con- 
clusions regarding the relative value of the different 
OAM series would be premature. 

6. Summary and Final Remarks 

We have examined more than 11 years of axial 
OAM signals calculated from output of a 1° model. 
Seasonal and shorter period signals are clear in both 
M r and Mq, but signals in Mq are generally stronger. 
Interannual variability is weak. There is a clear an- 
nual cycle, with amplitude nearly twice that of the 
semiannual cycle, and modulation from year to year. 
Atmospheric driving is found important for the pres- 
ence of subseasonal OAM signals. Strength of sea- 
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sonal signals varies considerably with region. Low lat- 
itudes contribute substantially to fluctuations in Afa, 
but there are other important signals over the South- 
ern Ocean and other regions. The ACC and related 
circulation in Southern Ocean are closely related to 
M r \ strong M r signals are also present at low latitudes 
but they are associated with recirculating currents 
at different latitudes and largely cancel out. OAM 
signals can explain some of the observed changes in 
LOD but their estimated amplitudes are substantially 
smaller than atmosphere-solid Earth residuals. 

It is interesting to note the similarity in the rela- 
tion between regional and global M r signals (banded 
structure in Figure 3b), and regional and global zonal 
wind stress torques (cf. Figure 5 in Ponte and Rosen 
[1993]). Previous works have shown the in-phase re- 
lation between seasonal signals in total M r and stress 
torque [PE94; B97]. Thus, both regionally and glob- 
ally, M r and wind torques vary nearly in phase. This 
behavior is rather surprising, since for angular mo- 
mentum balance a driving torque should lead OAM 
by 90°. Results of Ponte [1990] and B97 point to a 
plausible explanation. Seasonal wind stress and pres- 
sure torques nearly balance instantaneously, but the 
pressure torque lags slightly given the finite (albeit 
rapid) adjustment time in the oceans. Consequently, 
as can be easily seen by drawing a phasor diagram, 
the small resulting net torque should lead the wind 
torque by nearly 90°. (B97 noted this peculiar phase 
relation but left it unexplained.) Viscous torques fur- 
ther reduce the amplitude of this net torque but do 
not affect its phase (see Figure 12 of B97). From the 
relation between wind and net torques and from an- 
gular momentum conservation, it follows that M r and 
wind torque should be in phase as observed. The ar- 
guments have been presented for global quantities but 
should essentially hold regionally too [PR94]. 

From the comparisons between present and pre- 
vious results, it is clear that ocean models will give 
different axial OAM time series, depending on their 
formulations, forcing, bathymetry, etc. Using axial 
OAM series in the context of the planet’s angular 
momentum budget may provide a useful global con- 
sistency check on the models, much as proposed in the 
context of polar motion [Ponte et a/., 1998]. Compar- 
isons in a more controlled setting, as carried out, for 
example, by Hide et al [1997] for the atmosphere, 
should be most useful in this regard. 

Both Marcus et al [1998] and our results give weak 
OAM signals compared to the atmosphere-solid Earth 
residuals, at all frequency bands examined; similar 


findings are discussed by Ponte [1997] and Johnson 
et al [1999]. It remains to be determined whether 
this behavior is due to underestimating the strength 
of OAM contributions (xP and Xs~ A spectra are sim- 
ilar in shape and models are known to underestimate 
oceanic variability), to errors in X 3 and Xz (note 
the dropping with frequency of coherence amplitudes 
between x? and X 3 and the comparable power in 
Xz and X$ A series at fortnightly and shorter peri- 
ods), or to missing excitation sources not related to 
atmosphere or oceans. A more detailed examination 
of the axial angular momentum budget exploring all 
these issues is left for future work. In what regards 
improving OAM estimates and assessing their uncer- 
tainties, an important next step currently being pur- 
sued involves the use of state estimation methods to 
constrain the model to data. 
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Figure 1. Time series of Mq, M r , and M. For ease of comparison between curves, time means have been removed. 

Figure 2. Spectra and low-frequency periodograms of Mq (dotted lines), M r (dashed lines), and M (solid lines). 
Spectra are based on averaging over 20 adjacent frequency bands. First spectral point represents the 414-138 day 
band. Before Fourier transformation, series are tapered at both ends with a half-cosine bell. 

Figure 3a. Standard deviation of M r in kg m 2 s -1 . Contour interval is 10 22 kg m 2 s _1 . 

Figure 3b. Fractional covariance of M r with M r . Contour interval is 0.01. Light shading denotes negative values. 
Right panel shows fractional covariances added over longitude. 

Figure 4. Variability in ACC transport at Drake Passage calculated in OAM units as described in the text for 
(j> = 55° (solid line) and the globally integrated value of M r (dotted line). Small gaps in ACC transport series 
have been linearly interpolated. 

Figure 5a. As in Figure 3a but for Mq. 

Figure 5b. As in Figure 3b but for Mq. Contour interval is 2xl0 -3 . 

Figure 6a. As in Figure 3a but for M. 

Figure 6b. As in Figure 3b but for M. Contour interval is 2.5xl0 -3 . 

Figure 7. Time series of X 3 > X A > and X 3 (solid lines). Dotted line in bottom panel is the residual x?"'" 4 - Note 
the difference in scale. All time series have been detrended and high-pass filtered to supress variability at periods 
of 3.8 years and longer. 

Figure 8. (a) Spectra of X 3 (dotted-dashed line), X 3 (dashed line), xf A (solid line) and x 3 (dotted line) 
calculated as in Figure 2. (b) Coherence amplitude squared between X 3 and Xif+X? (solid line), xf and X3 (dashed 
line), and X3~ A and X3 (dotted line). Values above 0.15 are significantly different than zero at 95% confidence 
level, (c) Respective coherence phases in degrees. 
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Table 1 . Amplitude in 10 23 kg m 2 s -1 and Phase (Date of the First Maximum After January 1) of Mq, M r , and 
M for the Annual and Semiannual Periods. 


Annual Semiannual 



Amp 

Phase 

Amp 

Phase 

M n 

7.8 

Aug 27 

4.4 

Mar 15 

M r 

3.6 

Sept 1 

2.1 

Feb 24 

M 

11.4 

Aug 28 

6.2 

Mar 9 


x 1 0 kg 














x 1 0 kg 







CONTOUR FROM 1 TO 16 BY 1 






CONTOUR FROM 1 TO 1 6 BY 1 
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